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Abstract 

We examine the possibility of fermion mass generation in 2+1-dimensional 
gauge theory from the current algebra point of view. In our approach the crit- 
ical behavior is governed by the fluctuations of pions which are the Goldstone 
bosons for chiral symmetry breaking. Our analysis supports the existence of 
an upper critical Np and exhibits the explicit form of the gap equation as 
well as the form of the critical exponent for the inverse correlation length of 
the order parameter. 
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The possibility of dynamical generation of fermion masses is of funda- 
mental importance to our present view of quantum chromo dynamics (QCD) 
as a theory of the strong interactions. Central to our understanding of this 
phenomenon is the existence of a critical coupling. When fermions have a 
sufficiently strong, attractive interaction there is a pairing instability and the 
resulting condensate breaks some of the flavor symmetries, generates quark 
masses and represents chiral symmetry in the Nambu-Goldstone mode. This 
idea dates back to the earliest models for chiral symmetry breaking |l] and 
is prevalent in the modern literature 0. 

Recently the issue of critical coupling has been investigated in 2+1 di- 
mensional gauge theories || f|, ||. These theories provide toy models which 
exhibit simpler behavior than their 3+1-dimensional relatives || [7]] and are 
also of interest as effective field theories for some condensed matter systems 
[H]. Typically, their dimensionless expansion parameter is 1/Np, where Np 
is the number of quark flavors ||. Using Schwinger- Dyson equations in the 
1/Np approximation for QED and QCD, || have found that there is a crit- 
ical Np nt such that when Np < Np rit chiral symmetry is broken and when 
Np > Np rit it is not broken and quarks remain massless. 

In the case of QED, this result has been the subject of some debate 0, [|, 
H |lOj- Either an improved ladder approximation |5|, [10| or renormalization 
group computation |4| find no critical behavior and that chiral symmetry is 
broken for arbitrarily large Np. There are, however, numerical simulations 
[ PH of 2+1- dimensional QED which find an Np nt remarkably close to that 
obtained by ||. 

In this Letter we shall present further support for the existence of iV£ rit . 
We shall advocate a picture which is complementary to that of critical at- 
tractive quark-quark interactions and in which the critical behavior of the 
chiral symmetry breaking phase transition is governed by the fluctuations of 
the pions which are the Goldstone bosons for broken continuous flavor sym- 
metries. We shall argue that an upper critical Np is natural since the order 
parameter is renormalized by the N F /2 pions with classical coupling constant 
~ 1/NpNc where Nc is the number of quark colors. Their fluctuations are 
strong enough to destroy the ordered state when Np = Np Tlt ~ No- 

We are partially motivated by our recent study of strong coupling gauge 
theory |TJ| on the lattice. We showed that, in the strong coupling limit, a 
Hamiltonian lattice gauge theory with Nc colors and Np/2 lattice flavors 
of staggered fermions (because of fermion doubling this corresponds to Np 
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continuum flavors of 2-component spinors in 2+1-dimensions) is effectively a 
U(Np/2) quantum antiferromagnet, with representations determined by Nc 
and Np. We also identified chiral symmetry breaking with the formation of 
either commensurate U(l) charge density waves or SU(Np/2) spin density 
waves, i.e. Neel order. We found that the cases where Np/2 is an odd or 
even integer are quite different. 

When Np/2 is odd the strong coupling limit necessarily breaks chiral 
symmetry, no matter how large Np is and the condensate is a U(l) charge 
density wave. The ground state has a staggered structure with the SU (Np/2) 
representations given by the Young tableau with Nc columns and Np/4 + 1/2 
rows (and U(l) charge 1/2) on even sites and Nc columns and Np/4 — 
1/2 rows (and U(l) charge -1/2) on odd sites. (This was strictly true for 
U{1) gauge theory, and likely for U(N C ) and SU(N C ) gauge theory where 
it could be proved only with some additional assumptions about translation 
invariance.) 

We found that when Np/2 is an even integer, there is no U(l) charge 
density wave and the representation at each site was given by the Young 
tableaux with Nc columns and Np/4 rows. 

In either case of even or odd Np/2, there could be Neel order, which 
also breaks chiral symmetry. Quantum antiferromagnets with the kinds of 
representations we considered have been analyzed in [13j where they found 
that, for small enough Np, the ground state is ordered. Also, in the generic 
case, when Np is increased there is a phase transition with Np Tlt ~ Nc to 
a disordered state. In this picture, the large Nc limit is the classical limit 
where the Neel ground state is favorable and the small Nc and large Np limit 
is where fluctuations are large and disordered ground states are favored. 

For example, the SU(2) antiferromagnet with spin j corresponds to 4- 
flavor QCD with color group U(2j) and in particular to QED when j=l/2. 
It has a Neel ordered ground state for any j, corresponding to chiral symme- 
try breaking in the strong coupling limit of QCD. However, an SU(Np/2) 
antiferromagnet with Np a large multiple of 4 and in a representation of 
SU(Np/2) given by a Young tableau with a single column of Np/4 boxes 
corresponds to strong coupling QED with Np flavors of fermions. It is known 
that the ground state of this system is disordered with several competing flux 
and dimer phases || For some intermediate Np between 2 and oo the 
antiferromagnet has a phase transition where Neel order is lost. 

It is tempting to conclude that this critical behavior of antiferromagnets 
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is related to the critical behavior of continuum gauge theory found in P| 
using Schwinger- Dyson equations. In the following we wish to examine this 
question further in the continuum by analyzing the dynamics of the effective 
field theory for pions in the phase with broken chiral symmetry. 

Our lattice results seem to imply that there is a subtle difference between 
the cases where Np is an even or odd multiple of 2. It is not clear whether 
this is a fundamental difference in the continuum theory too, or merely an 
artifact of forcing gauge theory to live on a lattice. The numerical lattice 
simulations of QED in use staggered Euclidean fermions and can there- 



fore only study the case when Np is a multiple of 4. For compact continuum 
QED we argued in [12 that, if one uses Polyakov's idea |18| of taking the 



Georgi-Glashow model with spontaneous symmetry breaking SO(3) — > U(l) 
to obtain a photon with compact U(l) gauge group, then it is known |15[ 
that the minimal number of 2-component spinors the photons can couple to, 
consistent with gauge and parity invariance, is 4. Thus, if continuum QED is 
to be compact, Np is a multiple of 4. This could also apply to U(Nc) gauge 
theory because of the U(l) subgroup. However, the only such restriction on 
QCD with SU(Nc) gauge group constrains Np to be even. 

We shall consider QCD with gauge group SU(Nc) and Np (=2xinteger) 
flavors of massless 2-component quarks in Euclidean space 

s = Jd 3 xl—Y { f; v f; v + £ + A,)r\ (1) 

where e 2 is the dimensionless coupling constant and A is the ultraviolet cutoff. 
This action has U(Np) global flavor symmetry and also a Z 2 parity symmetry 
under the replacement (Ai, A 2 , A 3 )(x) — > (— A 1 , A 2 , A 3 )(x'), ip(x) — > 7i^(x'), 
— > —ip(x')'-fi with x' = (— Xi, %2i x%)- We shall assume that the ultravi- 
olet regularization preserves parity. 

An order parameter for the U(Np) x Z 2 symmetry breaking is the quark 
bilinear 

M ap {x)=$ a {x)ij) p {x) (2) 

M is a Hermitean matrix and transforms under U(Np) as M — > gMg^ and 
under parity as M — > —M. ( This is in contrast with its counterpart in 
3+1-dimensions, = ipiipR which is a complex matrix and transforms under 
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SUr(Nf) x SUl(Nf) as /i — > gixh). ) Flavor symmetry breaking is governed 
by the effective Landau-Ginsburg action 

S eff = J d 3 x tr (cxd^Md^M + c 2 M 2 + c z M A + . . .) (3) 

Note that, in the large Nc limit, the coefficients q in (|3]) are fermion loops 
with meson operator insertions which are naturally of order Nc [0] • 

We shall consder the symmetry breaking pattern U(Np) x Z 2 — > U(n) x 
U(Np — n). In this case, M has a constant vacuum expectation value 

M =< M >= const diag (1, 1,1,..., -1, -1, -1) (4) 

where there are n l's and Np-n -l's. We shall show that the case where n = 
Np/2, for which a residual parity symmetry can be defined, is dynamically 
favorable. 

We are interested in the dynamics of Goldstone bosons which are de- 
scribed by a sigma model with target space the Grassmannian 

U(N F ) x Z 2 



U{n) x U(N F - n) ' 
With the ansatz M(x) = g(x)Mog^(x) we obtain the sigma model 

S eff = / d3x ^jr tr ([gd,g\ M ] [gd^, M }) + S wz (5) 

where we have renamed the coefficient, which here plays the role of cou- 
pling constant and we have extracted its natural order in Nc- Swz is a 
Wess-Zumino term which must be added to the sigma model action in order 
to break an unwanted discrete symmetry and to obtain a sensible current 
algebra [0. (|) is equivalent to the gauged principal chiral model 

S cS = J d'x^f tr ((D 9 y ■ (Dg) + i\tfg - 1)) + S CS [V, W] (6) 

where D = d + i(V + W) and V and W are Hermitian gauge fields which 
have components in the upper left n x n block and in the lower right (iVV — 
n) x (Np — n) block respectively. In (|6]) the Chern-Simons action is 

S CS [V, W}=i— [ d 3 xti (vdV - WdW + -V 3 - -W 3 ) (7) 
4:7i J V 3 3/ 
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and we have introduced an Np x Np Hermitean Lagrange multiplier field A 
to enforce the constraint gg> = 1. Eliminating V and W in using their 
equations of motion yields (|5]), up to higher derivative terms which come from 
approximating the Wess-Zumino term by the Chern-Simons action. In |17| it 
was argued that, the model described by (@) has features remarkably similar 
to the commonly accepted features of the Skyrme model of 3+1-dimensional 
QCD, such as solitons which have Fermi statistics and behave like baryons. 
They also argued that to obtain the correct current- current commutation 
relation, one must set 9 = Nc- To be general, we shall keep 9 arbitrary in 
the following analysis. 

In order to study the quantum properties of this model, we first note that 
the field g appears quadratically and can be integrated to get the effective 
theory 

S eS = N F TR \n(—D 2 +iX)-J d z x % -^- tr A + S CS [V, W\ + S ghost (8) 

where, to fix the gauge freedom and properly define the quantum problem 
we have added the Faddeev-Popov ghost action, 



5 gh0 st = / d 3 x tr (j^(dV) 2 + jp(dW) 2 + dc\d + iV)c + 



dd ] {d + iW)d 



As in the standard approach to sigma models [19| , the remaining 
integral over A, U and V is done by saddle-point approximation. It is assumed 
that the gauge fields are zero at the saddle point. The saddle point value of 
A provides a mass for the chiral field g in @. When this mass is non-zero g 
fluctuates about g=0 and the model is disordered. When the saddle point is 
at zero, we obtain the ordered phase. Assuming a constant saddle point and 
putting iX = fi 2 , we obtain the gap equation 



f 2 2vr 2 I A /i 



A solution of this equation exists and the model is disordered if 

N F > 2n 2 ^ (10) 
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with the equality giving the condition for criticality. For fixed f 2 and Nc 
we can interpret this as an equation for critical Np. When Np exceeds 
N F nt = 2ir 2 Nc/f 2 , the chiral symmetry breaking condensate is unstable. 

Note that, unlike the case of the more conventional O(N) non-linear sigma 
model, the saddle point approximation in the present case is not controlled by 
any small parameter such as 1/N. This is because, even though the effective 
coupling constant of the saddle-point method is ~ 1/iVp, there remain of 
~ Np degrees of freedom. The correct large Np expansion of (^) would 
involve the topological expansion where, like the large Nc expansion of QCD 



| fL6f| , it is necessary to sum all planar graphs to obtain the leading order. 

In spite of this, we have two reasons to believe in the validity of fllQf) . 
First of all, as we shall see in the following, corrections from the next order 
in the saddle point approximation are indeed smaller than the leading order. 
Secondly, alternative to the saddle point analysis, we could do a weak cou- 
pling expansion. If Np is large enough, criticality in ( |T0D is obtained in the 
weak coupling region. 

The computation of the gap equation to next order is somewhat sophisti- 
cated. We must expand (M) to quadratic order in V, W, A, drop linear terms 



2C] and perform the functional integral in the Gaussian approximation. 



at3t , f A d 3 k , .,, 2 , ANcNpfi 2 2 -r j- f A d 3 k 1 



where 



+N F J tr (-AAA + -F V UF V + -F W UF W + ...)+ S cs + S ghost (U] 

\k\ 




— : — r arctan ■ . 
47r|A;| 2\/j, 

k 2 + Aii 2 
H 7T7Z — arctan 




k 2 \k\ 3 \2\fi\ 

V is the volume, we have set A — > — ifi 2 — i5fi 2 + A and 5fi 2 = —2A 2 /tc 2 is a 
counterterm which arises from an additive shift of A and which is necessary to 
cancel a quadratic divergence. This infinite shift of the Lagrange multiplier 
field is a standard feature of sigma model renormalization |TJ| . The result of 
integrating the Gaussian fluctuations is 
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Nl , A n 2 + (N F - n) 2 , /i9tt9 ,„ . , 9 \ 
+ In A + J - In [k 2 U 2 + {d/Aix) 2 ) + 

1 2 l + 2n 2 /N 2 -2n/N F ' 




N 2 Afi 2 - 



19/ir 



2\2 



- I . A 16 1 + 2n 2 /A^ 2 - 2ra/A^ F A 



In — + /„, NO - In— + (12) 

3 V vr /ia 7T 4 l + (8#/vr) 2 ^ v ; 

Here, a and /? are (unknown) constants. The method for computing the 
integrals in described in the Appendix. If interpreted as an effective potential 
for ji, S C R is upside-down and apparently unstable. This originates with 
subtleties in dealing with complex saddle points. We refer the reader to the 
standard literature on the subject pl |. The corrected formula for the critical 



line is 

^crit ,j^ c / 2 4(l + 2n 2 /iV 2 -2n/iV F ) \- 1 

Nf =27C l 2 -{ ~ 2 + 1 + J 

When 9 (which is equal to Nc) is large, the right hand side differs by about 
20 percent from the leading order estimate. The largest Np nt occurs when 
n = Np/2. Thus, the first and therefore most stable ordering which occurs 
as as we lower Np is the parity symmetric phase n — Np/2. 

We can also interpret the logarithms in (|T^) as changing the exponent, 
| /j, | 3 — ► |yu| 3< - 1+7 \ where 

8 / TT 2 l + 2n 2 /N 2 F -2n/N F \ 
7 ~ 3^ y 1 ~ T + 1 + (80/tt) 2 J 

and the gap equation has solution 

/J, ~ (Np — ATcrit)V(2+3 7 ) 

which exhibits the critical exponent (again a small correction of the leading 
order result 1/2) for the inverse correlation length of the order parameter. 

In conclusion, we note that the beta function of the Grassmannian sigma 
model has been computed using the e-expansion about 2 spacetime dimen- 
sions in ||22|| . In their notation it is 

P(t) = (d-2)t- N F t 2 - ... 
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and has an ultraviolet stable fixed point at t = {d— 2)/Np- In our notation 
the coupling constant is t = f 2 /27i 2 Nc, and in the leading order, the critical 
point ( |T(1| ) occurs at the zero of the beta-function. This is an ultraviolet stable 
fixed point (it cannot be infrared stable since the flow to low momentum 
should be interrupted by mass generation). 
Appendix: Calculation of Integrals 

As an example of the integrations needed to find the quantum corrections 
to the gap equation, consider the integral J[A//i] = J^ 2fl dx x 2 In- arctan x 
which we want to evaluate in the limit A/fi —>■ oo. To evaluate the integral, 
we first extract the divergent parts as 

f A/2 ^ , 2 2 2 1 2tt 2 -8 \ 

I A/// = / ail In- arctan x H h - „ „ 9 , r + 

jo \ 7r ttx ir z x z 6ir' > x z {x + a) J 

A 2 A 2tt 2 - 8 , . . /n . 

+ ^— ln(A/2/xa) 



47T/X 2 7T 2 /i 37T 3 

In the first term the infinite cutoff limit can safely be taken, to get a function 
of a. However, a is an arbitrary number, so the integral can be parameterized 
as 

rri , , A 2 2A 2vr 2 -8 , .. , n 
WfA = ~~a 2 " — + In A/2//a 

where a is a fixed, but unknown (and irrelevant) constant. 
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